function [coefficients, fitted_curve] = polynomial_fit(x, y, order)
    % 多项式拟合函数
    % 参数：
    %   x: 自变量数据点
    %   y: 因变量数据点
    %   order: 拟合多项式的阶数
    % 返回值：
    %   coefficients: 拟合多项式的系数
    %   fitted_curve: 拟合曲线
    
    % 使用最小二乘法进行多项式拟合
    coefficients = polyfit(x, y, order);
    
    % 根据拟合得到的系数创建拟合的多项式函数
    fitted_curve = polyval(coefficients, x);
    

end
